Parasitic pumping currents in an interacting quantum dot 
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We analyze the charge and spin pumping in an interacting dot within the almost adiabatic limit. 
By using a non-equilibrium Green's function technique within the time-dependent slave boson ap- 
proximation, we analyze the pumped current in terms of the dynamical constraints in the infinite-?/ 
regime. The results show the presence of parasitic pumping currents due to the additional phases 
of the constraints. The behavior of the pumped current through the quantum dot is illustrated in 
the spin-insensitive and in the spin-sensitive case relevant for spintronics applications. 
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I. INTRODUCTION 

The quantum pumping proposed by Thouless- 1 is a 
phase coherence effect able to pump dc current by the 
out-of-phase adiabatic modulation of (at least) two sys- 
tem parameters in absence of a voltage bias. The phase 
difference ip of the external signals produces a charge 
current proportional to wsin((^), uj being the adiabatic 
pumping frequency. After the Thouless theory, Brouwer- 
formulated a scattering approach to the adiabatic quan- 
tum pumping showing that the dc charge current pro- 
duced by the pump is related to the parametric deriva- 
tives of the scattering matrix of the device. 

Following the Brouwer formulation several studies, 
both of experimental^ and theoretical nature, have been 
performed. The theory of quantum pumping for nonin- 
teracting systems^ has been developed in several direc- 
tions. In case of interacting electrons the computation 
of the pumped charge becomes rather involved and few 
works have addressed this issue for different systems^ and 
in specific regimes. As for the case of interacting quan- 
tum dots, the pumped charge in a period was calculated 
by Aono^ by exploiting the zero-temperature mapping of 
the Kondo problem within non-equilibrium Green's func- 
tion (NEGF) technique. A very general formalism was 
developed in Ref . Q based on the adiabatic expansion of 
the self-energy in the average-time approximation, while 
a linear response scheme was employed in Rcf.|8]. An- 
other interesting study 9 was performed aiming at gen- 
eralizing Brouwer's formula for interacting systems to 
include inelastic scattering events. Alternatively to the 
NEGF approach, a powerful diagrammatic technique to 
treat the interacting quantum pumping is constituted by 
a generalized-master equation method^. 

Very recently the idea that the presence of an hidden 
classical variable inside a scattering region can strongly 
affect the dependence of the current pumped through the 
system has been proposed. This effect is produced by 
the dynamics of the classical variable which modifies the 
standard phase relation introducing a dynamical phase 
shif t 11 ' 12 accompanied by parasitic contributions to the 
pumping current. 

In this paper we develop a time dependent mean field 



slave boson theory for an interacting quantum dot follow- 
ing a non-equilibrium formulation a la Keldysh. In par- 
ticular, we derive a time-dependent mean field equation 
for the slave boson and the constraints which allows us 
to study the strongly interacting regime of the pumping 
beyond the adiabatic (i.e. zero frequency) limit. How- 
ever, since the cost of the analysis is growing with the 
degree of non-adiabaticity, in this work we focus on the 
almost adiabatic limit (i.e. pumping frequency small but 
finite). In particular we focus on the effect of the dy- 
namical phase of the slave boson on the pumped current, 
evidencing features not discussed before. Furthermore, 
the proposal of a quantum pump able to produce spin 
current without polarized electrodes is also discussed. 

The organization of the paper is the following: In 
Sec|nl we introduce the model Hamiltonian and derive 
the free-fermions-like effective Hamiltonian by means of 
the slave boson treatment. In Sec lIIII we derive the rele- 
vant Green's functions within the Keldysh approach and 
a formula for the pumping current is derived within the 
single photon approximation. In Scc lIVI we analyze the 
constraint equations and propose an approximate solu- 
tion suitable within the almost adiabatic limit. In SecfVl 
we discuss the results obtained both in the case of a spin 
independent quantum pump and in the case of a spin 
sensitive device able to generate spin current. The con- 
clusions are given in Sec lVIl 



II. 



THE MODEL AND FORMALISM 



We consider a quantum dot (QD)-based pump cou- 
pled to noninteracting leads. A particle flux is generated 
through the QD by means of an almost-adiabatic mod- 
ulation of two out-of-phase gate voltages controlling the 
transparency of the tunneling barriers. 
The Hamiltonian of the quantum pump depicted in Figfj] 
can be written as follows: 
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where the operator c kaa creates an electron of momen- 
tum k and spin a in the lead a, while S is the creation 
operator of an electron state on the interacting quan- 
tum dot, e being the energy of a single occupied electron 
state. The electron-electron (e-e) interaction energy U, 
described by the third term, has to be considered when 
the dot screening length X s is bigger than the typical size 
of the quantum dot. The last term describes the tun- 
neling between the quantum dot and the leads. In this 
work we are interested in describing the strongly inter- 
acting limit (U — > oo)s*. In the intermediate {/-regime 
the Hilbert space of the isolated QD has electron states 
{|0}, |tr), |2}} (empty, single and doubly occupied) which 
are created by application of the creation operator to the 
empty state |0), i.e. \a) = 4|0>, |2> = | U) = d\d\\0) 
(notice that the order is arbitrary but have to be fixed 
since | |4) = — | it))- I n the previous basis the relevant 
operators can be written asi^: 

dl = |a')(0|+a'|2)(a'| (2) 

d a = \0)(a\+a\a)(2\ 

n a = |a}<a| + |2)(2| 

n tni = |2}(2|, 

and the canonical anti-commutation relations 4} = 
&iiv> {dn,d v } — 0, {4>4} = can be obtained. In case 
of strong interaction (U — > oo) the occupation of the 
|2)(2| subspace becomes strongly suppressed and thus 
the commutation relations^ become {4; 4) = 5^ u — 
yLvdXdp,, where fi,v € {til} = 1} an d the effec- 

tive completeness relation is 1 = |0)(0|+^ cr \<j)(a\. Thus 
the electron operators, with the exclusion of the doubly 
occupied state, take the form d a = |0)(rr|, 4 = |cr)(0|. 
By introducing the reference state \ref) and the opera- 
tors b and f a we define |0) = tf\ref) and \a) = fl\ref) 
and thus the creation/annihilation operators can be writ- 
ten in terms of bosonic (i.e. b) and quasi-fermionic (i.e. 
f a ) operators^ as: d a -> tff a , 4 -» /+6. Within the 
slave boson representation the Hamiltonian is: 

H SB = H(d a 6t /a . d t _> } \ h) + \(tfb + J2 flU - I)- 

(X 

(3) 

where the single occupation constraint is included by the 
Lagrange multiplier A which is fixed by the equation: 

dx(H SB ) = -> + <6 + 6> - 1 = 0, (4) 

while the slave boson operator b evolves in time according 
to the equation of the motion ihd t b = [b, Hsb], thusi£: 

ihd t b = \b + Y J V ka * (t)4 aa J (j . (5) 

kacn 

We treat the slave boson operator b within the mean 
field approximation ((b) = B, (w) = £>*) and the original 
problem containing strong correlations is replaced by a 



constrained free-fermions-like theory whose dynamics is 
completely described by the following Hamiltonian: 

Hsbmf = H leads + J2(t + Ht))ftf°+ (6) 

a 

kaa. 

+ \{t){\B{t)\ 2 -i), 

with the dynamical constraints: 

ihd t B(t) = \{t)B{t) + VkUtMaaU) (?) 
kaa 

J2(rtU) + Mt)\ 2 -1 = 0. (8) 

As shown in (6) the interaction produces a renormal- 
ization of the QD energy level (e — > e + A(i)) and 
a redefinition of the tunneling amplitudes (Vfc acr (i) — > 
Vfcatx (£)£>*(£)). This aspect is very important in the con- 
text of the quantum pumping. In fact, since the time evo- 
lution of the slave boson field is governed by the equation 
([7]) , a dynamical phase shift with respect to the phase of 
the external driving signals is expected. Within the adia- 
batic case, all the phases are rigidly related to the exter- 
nal driving signals, while at finite frequency this relation 
is lost and interesting contribution of hidden phases can 
emerge. 



III. NEGF METHOD 

In this section we employ the non-equilibrium Green's 
functions (NEGF) formalism 1 — to derive the current flow- 
ing through the system when subject to an almost adia- 
batic quantum pumping (small but finite pumping fre- 
quency uj) in which we modulate the tunneling rates 
r?(t) = r« + r? j(1> sin(wt + y a ) defined by T«(i) = 
27rp a \Vk F aa(t)\ 2 , where p a is the density of states at the 
Fermi level kp. 

Due to the free-fermions like form of the Hamiltonian the 
retarded GF of the quantum dot (QD) uncoupled from 
the external leads is (h = 1): 

g r sp {t,t / ) = -i5 sp 9{t-t')exp{-i J dh(e + X(h))}, (9) 

being s,p E the spin index. When the QD is cou- 

pled to the leads, we must take into account the transition 
rate of an electron through the system. The quantity en- 
coding this information is the retarded self-energy which 
can be written within the wide band limit (WBL) as: 

K P (ti,h) = -iS sp 6(h-h) ]T -^pl\B(h)\ 2 , (io) 

a— l,r 
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and thus the retarded Green's function of the QD coupled 
to the leads takes the form: 



G r s Jt,t')=g r 8 Jt,t')eM 



r dtl ^ mA mtl)n 

Jv a=l,r 



(11) 

which depends on the renormalized quantities |£>(t)| 2 and 
A(t). The corresponding advanced quantities, i.e. G a and 
E a , are computed directly by using the general relation 
S",(ti,ta) = ^a(*2 S ti)*, S being G or E. Finally, the 
Langreth rules^ can be employed to compute the lesser 
self-energy as 



E< {h,t 2 ) = iS sp f(ti - h) 



r«( tl )\B( tl )\ 2 . (12) 

a— Z,r 

It is proportional to the Fourier transform f[t\— ta) of the 
Fermi function /(2£) (notice that the chemical potential 
is the same for both the leads). 

In order to calculate the current flowing through the 
QD we need to calculate the lesser GF Gf p (t\ : t 2 ) of the 
QD exploiting the Keldysh equation^ and by using G r l a 
and E < . Together with the single particle Green's func- 
tion we need yet to solve the constraint equation. Ex- 
plicitly we can rewrite Eq.([7J and <j8j) as: 



id t B(t) = [\(t) + Y,G< a (t,t) 



E 



r„(t) 



B(t) 



(13) 



\B{t)\ 2 = l+»^G<(t,i) (14) 

Let us note that the interaction effects are taken into 
account by the constraints. This aspect becomes evi- 
dent considering the limiting case in which = T (spin 
insensitive tunneling). Solving (fT4")l for ^2 G^ a and sub- 
stituting into (jT3"]l . we obtain 

.i^±\B(t)\ 2 B(t) + ..., (15) 



id t B{t) 



[X (t) + l T M]B(t)- 



in which the nonlinear term \B\ 2 B accounting for the in- 
teracting nature of the problem appears. The obtained 
equations are similar the one derived by Langreth et al. 
in Ref. [22j even though in our time-dependent theory we 
expect an easier manipulation. 



A. Non-equilibrium currents 

Within the non-equilibrium theory the spin resolved 
current pumped in the lead a can be written in terms of 
the QD Green's function in the following form: 



2e 



Re 



dE 1 dE 2 dE 3 



h (2tt) 3 
x eMKE3-E 1 )t][G r aa ,(E 1 ,E 2 )^< a , a (E 2 ,E 3 ) 

+ G<A&, &)]}> 



(16) 



where the two-times Fourier transform^ of the GF have 
been introduced. In this work we are interested in de- 
riving the pumping current generated by an almost- 
adiabatic perturbation and thus we expand all the rele- 
vant time dependent quantities by considering only one- 
photon contributions as follows: 



A(t) 

B(t) 
ZUE l7 E 2 ) 



71=0, ±1 

E 

n=0,±l 



X n exp{— inuit} (17) 
r ";nexp{-ina;i} 
B n exp{— inwt} 

n=Q,±l 

S sp Zi, v (Ei)6(E 1 -E 2+V u), 

J1=0,±l 



where S is G or E, while £ £ {r, a, <}. Within the 
one-photon approximation, the spin resolved dc cur- 
rent J£ pumped in the lead a takes the form 1° = 
(2e/h) J2 v =±i ^ e [-^Ir:»)] where we defined the following 
quantity: 



^[G r a -ri{El)^'aa--ri 



(2tt) 3 



(E 1 + 7]L0) + (18) 



G< v {E l )T. a aa ,_ ri {E l + V u)], 



which contains information about the pumping cycle and 
the absortion/emission processes of one photon. Within 
the single photon approximation the behavior of the 
pumped current with respect to the QD energy level is 
determined by the two integrals: 



1^ = jJ^f{E)V°{E)V°{E) 
K n = JJ^f{E)V%{E)V^EY, 



(19) 



where 



V°{E) = [E + r)u;-(e+ A ) + i{ la /2)\ 



(20) 



while Ao represents the static part of the Lagrange mul- 
tiplier and 7o- is a renormalized linewidth depending on 
the slave boson field: 



E 



[\B Q \ 2 Y% + 2Re{\B \ 2 K :1 + (21) 
+ TZMBt+BoBlJ} . 



la 



As shown in (|2Tj) the finite frequency effects modifies the 
mean lifetime of the electron on the QD due to the in- 
terference effects originated by the slave boson field B±\. 
Furthermore, since picks up a phase due to the ex- 

ternal driving signals T". ±1 — ii-^ 1 exp{T-i<^a}> the GF 
of the problem depends separately on ip a and not via the 
phase difference tp = tpn — tpt,. The last observation im- 
plies that the quantum state of the dot is affected by the 
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and 



Left lead 



FIG. 1: The system described in the main text: A Quantum 
dot (dashed region) coupled to external leads via tunnel bar- 
riers whose transparency can be modulated in time via top 
gates. 



precise external phases which couple the system and not 
only by the phase difference tp applied to the nanostruc- 
ture. Moreover, since the slave boson field is a complex 
quantity, S±i,o = |£>±i,o| exp{i0^' }, interference terms, 
i.e. of the form cos(</>^ — (j>°) appears in T>°{E). Also 
additional phases associated to the constraint evolution 
appear. Thus beyond the almost adiabatic limit, the ex- 
istence of a great number of dynamical phases (and the 
corresponding interference terms) provides a decoherence 
source for the electron transport. 



IV. CONSTRAINTS AND THEIR 
APPROXIMATE SOLUTION 



Constraints within the single-photon 
approximation 



(23) 



where the auxiliary functions <&( m > n )(E) involving the 
GFs of the QD are reported in the appendix [A] In prin- 
ciple the finite frequency problem gives us an infinite set 
of coupled constraints, while within the one photon ap- 
proximation we obtain a reduced number of equations. 
However, from the numerical point of view, the solution 
of Eqs. (l2"2"|) -(f2"3" f may be in principle a non-trivial task 
since there can be numerical instabilities in the determi- 
nation of the solution. When an almost adiabatic regime 
is considered (retaining a finite frequency u) one consid- 
ers the static constraints as perturbed by the remaining 
ones which on their turn can be treated by using some 
suitable approximations to be checked at the end of the 
computation. 



B. Approximate solution of the constraints 

In this subsection an approximate solution of the con- 
straints shown in Eqs. (|2"2"|) - (|2"3"|) is provided within the 
almost adiabatic case. From the analysis of the paramet- 
ric expression of the current pumped through the system 
one can observe that the relevant quantities to be com- 
puted are: A , Ai, \B \ 2 and BqBi + BqB*_ x . Since we 
are modulating the tunneling amplitudes via top gates 
the contribution to the current coming from the modu- 
lation of X(t) can be considered as a parasitic pumping 
effect whose intensity is very weak when Ai — > 0. Thus 
we need only to produce an approximate expression of 
A±i and B^B\ + B B*L 1 to the first order in r ff; i. The 
constraints which need strictly to be valid are: 



The expression of the current pumped through the sys- 
tem can be calculated once the GF is determined. Since 
it depends parametrically on {Ao, A±i; Bo, B±i} one has 
first to solve the constraint equations. By exploiting 
the Fourier expansions of Eqs. p^| - p^|) within the one- 
photon approximation we get the following set of con- 
straints: 

A £ = -$ (0;0) (£ = 0) (22) 

(Ao - u)Bi + XiB = - * (1_W ~ / * ) ( J S = w) 
m=o,i 

(Ao+w)B_i + A_iBb = - ^> { ~^^(E = -u), 



— f{E)Re{V°{E)} (24) 
|£ | 2 = l-]>>c)| 2r ^ / ^f{E)\V°{E)\\ 

The previous equations come from the first equation of 
each set of constraints given in Eqs. (l2"2"l) - (|2"3"|) . They have 
the same structure as in the adiabatic case£ (not pre- 
sented here), even though, differently from that case, an 
hidden dependence on B±\ is present in T>q(E) and thus 
higher order constraints have to be considered in order 
to solve the problem. In particular, by considering the 
second equation in (f2"3")l it is possible to formally solve it 



5 



in the following form: 

B^Bi + B B*_ i = (25) 

i + E.r^o/ffwmW^) ' 

where the right hand side of the previous equation de- 
pends on Ao, |So| 2 and B±i, while we defined f^'(E) = 
f(E±Lj). Expanding the r.h.s. of Eq. (j2"5")) up to the first 
order in r a -x and taking the limit w — »• 0, the dependence 
on B±i disappears and thus Eqs. (|24)) - ([25|) can be solved 
independently from the other equations. Following a sim- 
ilar procedure, considering the limit r CT; o/Ao -C 1 we ob- 
tain: 



Ai 



Ao \ -< 
2 ^ 



B \ 2 r, 



<x:l 



dE 



V%(E)\ 2 f(E) 



(26) 



V. RESULTS 

A. Spin-insensitive quantum pumping 

In the following analysis we measure the energy in unit 
of To = 10 fieV, while the current is normalized to the 
quantity io = — eTo/h (io « 2.56 nA). The zero tem- 
perature limit is assumed. In Figj2]we show the current 
pumped in the right lead Ir as a function of the position 
of the bare level of the dot e and by setting the remaining 
parameters as follows: T^ L ^ R = 1.5, T^u" = 0.05, 
lu = 0.05, ipL = 0, ifiR = 7r/4. Since we are in the spin- 
insensitive case each spin channel contributes the same 
quantity to the charge current Ir (Ir = Ir)- The lower 
curve in Fig[2] represents the contribution to the current 
induced by the modulation of the renormalization of the 
energy level on the QD (i.e. A(t)), the upper curve is the 
current pumped by the modulation of the tunnel barriers, 
while the middle dashed curve represents the total cur- 
rent which is sum of the previous terms. When e>/i the 
current induced by X(t) is a vanishing quantity and the 
total current is dominated by the modulation of the tun- 
nel amplitudes; on the other hand, when the energy level 
of the QD e is placed below the Fermi energy fj,, the main 
contribution to the current is given by the modulation of 
X(t). From the analysis above we observe a competition 
between the two different pumping mechanisms which in- 
duces a change of sign of the total current and presents 
maximum values in the pA range. The behavior of the 
current shown in Fig[2] can be understood by analyzing 
the slave boson parameters as a function of (e — fj,)/To. 
Indeed in FigJ3]we plot the £>o| 2 vs (e— //)/Fo curve fixing 
the other parameters as in FigJ5] As expected the elec- 
tron density on the QD, roughly proportional to 1 — £>o| 2 , 
becomes a vanishing quantity when the energy level e is 
well above the Fermi energy. In this situation a weak 
renormalization of the tunneling amplitudes is observed, 
while in presence of a non- vanishing electron density on 
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FIG. 2: Total Charge current Ir (dashed line) as a function 
of (e — aO/To computed setting the remaining parameters as 



follows: r 



a=L/R 



1.5, r: 



t=L/R 



0.05, u = 0.05, <p L = 0, 



ifR = 7r/4. The lower curve represent the current induced by 
the modulation of A, while the upper curve depends only on 
the modulation of the tunneling amplitudes. 
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FIG. 3: Holes density |So| 2 as a function of (e — aO/To com- 
puted setting the remaining parameters as in Fig[2] 



the QD the tunneling processes are strongly renormal- 
ized. A similar effect is evident in Figf4]where we present 
Ao as a function of (e— /z)/To. As can be seen the effect of 
the constraints becomes more pronounced at decreasing 
(e — fi), i.e. when approaching the single occupied state 
of the dot. 

Another effect of the renormalization of Ao is shown 
in FigfS] where the absolute value of the oscillation am- 
plitude of A(i) is plotted as a function of the QD bare 
level. As shown in that figure, below the Fermi energy 
(i.e. close to e — fi w — 7Fq) |Ai| becomes comparable 




FIG. 4: Ao as a function of (e — h)/Tq computed setting the 
remaining parameters as in Fig[2] 



6 



to the oscillation amplitudes of the tunneling rates T" 
(which approach the value ~ 0.05r ) and thus the pump 
starts to feel the effect of the modulation driven by the 
internal dynamics of the system. Such effect is evident 
if we compare FigfJ] and FigJSJ Indeed FigfJ] shows that 
the A-induced current (lower curve) is enhanced when 
Ai| (see FigJS]) is close to its maximum value within the 
parameters range considered. In Fig. [5] we report the 




FIG. 5: | Ai| as a function of (e — fi)/To computed setting the 
remaining parameters as in Fig[3] 
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FIG. 6: Charge current II (upper panel) and Ir (lower panel) 
as a function of the pumping phase ipR computed fixing the 



remaining parameters as follows: T 



a = L/R 
<t;0 



i.5, r 



a = L/R 



0.05, u) = 0.05, <pl = 0. For both the upper and lower panel 
the dashed curve is computed taking (e — /-i)/^o = 1, the 
dashed-dotted curve for (e — f^)/^o = — 1, while the full-line 
is adopted for (e — fi)/To = 0. 



dependence of the charge current on the pumping phase 
tySfl . It exhibits a complicate behavior due to the presence 
of phase differences, namely (j)®'^ and </)^, induced by the 
slave boson field and by the modulation of A(i). Due to 
the presence of these phases, the external phase differ- 
ence <pr — ifL is not the only relevant phase, differently 
from the noninteracting case. Indeed, in the infinite-C/ 
limit terms of the form sm(ip a ± ((fy^ — 0°)) appear in 
the charge current. The interaction-induced terms are in 
general of the same order of magnitude of the ones usually 
obtained in the free-electron quantum pumping and thus 
strongly affect the usual sin(^)-behavior. Although the 
charge current contains terms having the usual forrn^ 
C+Asm(tpB — ipl)+B cos(ifR— ipl), its general aspect for 
the parameters range considered here is dominated by a 
cosine-behavior typical of a rectification contribution, i.e. 
of the form A + B cos((fu). The above result is very sim- 
ilar to the one found in Ref.[25| where the nonadiabatic 
noninteracting case has been considered. In the almost 
adiabatic case treated here, the external ac signals drive 
the system under non-equilibrium condition inducing an 
oscillation of the internal fields B(t) and X(t). The latter 
situation reproduces the external forcing pattern similar 
to the one proposed in Rcf.[25J. The Il/r vs ipa curves 
present a cosine-dependence, i.e. a rectification contri- 
bution, while the precise shape of the curves depends on 
the energy e of the QD (see the dashed, dot-dashed and 
full line in FigjBJ). To render the analysis more quanti- 
tative, while maintaining the structure of the equations 
simple, we explicitly report the expression of the current 
pumped in the lead a omitting the current term related 



to the modulation of A: 

i: - -(e/h) ]T \^[Re{R^}{r a r a , + i a i a ,)\27) 

r}=± a' 

+T]Im{K2 n }(r a L, s - r a i Q ) j + (...), 

where a — L8 a R + R^aL, while the spin-dependent coef- 
ficients r a and i a are defined as follows: 

r a = sm(y a ) rg; "f o1 +rg. |Bbl|gi| £ cos(0g - 0^8) 

la = cos( ya ) rg; "f° |2 +r° |g ||gi| £ vM$ - 4)- 
1 v =± 

As explicitly shown, the charge current contains all the 
phase differences (i.e. internal or external phase differ- 
ences) inside the scattering region. Thus in the strongly 
interacting case the presence of new phase differences is 
the main modification compared to the free-fermion the- 
ory. 



B. Spin-selective quantum pumping 

In this section we study the quantum pumping through 
the system shown in Fig [7] It consists of a strongly inter- 
acting QD connected to two non-magnetic leads via thin 
magnetic tunnel barriers. The magnetic barriers work 
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FIG. 7: The spin-sensitive quantum pump described in the 
main text obtained by connecting via spin-filter barriers (SF) 
two non-magnetic leads with an interacting quantum dot 
(QD). An ac modulation of the barrier strength is allowed 
via top gates Gf/G2. 



as a spin-filter (SF) whose effect is controlled by the po- 
larization p a £ [0, 1] of the a-th tunnel barrier. Thus, 
the spin and charge currents generated by the quan- 
tum pumping procedure are strongly affected by relative 
orientation (parallel or anti-parallel) of the SF barriers. 
Similar to Refs . (26|j27| . we take the tunneling amplitudes 
of the following form: T^(t) = T a (t)(l + p a a), where 
the polarization p a is assumed to be a time-independent 
quantity, while the tunneling amplitude T a (t) is taken: 



v a (t) = 



exp(^ a ) 
2i 



exp(iwt) 



c.c. 



(29) 



In the presence of the ac modulations of the gates Gl/ G2 
and thanks to the spin-selective properties of the SF, 
a spin current I s = al a is generated other than a 
charge current I c = ^2 a Ia- The proposed mechanism 
which produces a spin polarized current by only using the 
filtering properties of the magnetic barriers and avoid- 
ing the injection of spin polarized electrons by ferromag- 
netic leads, is very similar to the one proposed in Ref.[28( 
within the dc case and is conceptually new. In FigJ5]we 
study the charge (full line) and spin (dashed-dotted line) 
currents as a function of the bare energy level e of the QD 
setting the remaining parameters as in the figure caption. 
For non-vanishing polarizations (i.e. p L = pa = 0.4) of 
the SF barriers a sizable spin current is observed. Fur- 
thermore when e ~ fi — 5Tq the charge current takes neg- 
ligible values and the pump works as a pure spin current 
injector. A relevant observation is that the pure spin cur- 
rent obtained here is not related to the Zeeman energy of 
the QD (not present here) but is originated by the spin 
filtering properties of the magnetic barriers connecting 
the QD to the external unpolarized leads. Since the only 
relevant parameters are the polarizations Pl/r induced 
by the SF barriers, in Fig|9]we study the charge current 
as a function of the relative orientations of pl and pa by 
setting the remaining parameters as follows: Tq =l = 1.5, 

T«= R = 1.5, TZ =L/R = 0.3, u = 0.05, ip L =0,tp R = tt/4, 
(e — (J>)/Tq = 0. The analysis of FigJH] shows the charge 
current vs the polarization of the SF barriers. It presents 
high values of the current within the parallel configura- 
tion, while a suppression of the current is observed for 
the anti-parallel case. This behavior is explained by ob- 
serving that for parallel polarizations of the SF barriers 
only one spin channel is perturbed by the scattering off 
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h* 0.001 
0.000 
-0.001 
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FIG. 8: Charge current I c (full line) and spin current I s 
(dashed-dotted line) pumped in the right lead as a function 
of (e — fij/To by setting the remaining parameters as follows: 

v%= L = i.5, r%= R = i.3, rz =L/R = 0.3, w = 0.05, <p L = 0, 

fa = 7r/4, pl = Pr = 0.4. 
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FIG. 9: Contourplot of Charge current I c as a function of 
the polarizations Pl/r by setting the remaining parameters 
as follows: rg= L = 1.5, T^ R 
<PL = 0, ip R = 7r/4, (e - n)/T = 0. 



1.5, rZ~ L/R = 0.3, u = 0.05, 



the magnetic barrier, while for the antiparallel case both 
the spin channel are perturbed leading to a less intense 
current. The analysis of the spin current I s performed 
in FigJTU] shows that the relation I s w Ap L + Bp^ with 
A ?» 0.94 x 10~ 3 and B 3.3 x 10~ 3 is obeyed for the 
considered parameters. The latter observation implies 
that an intense spin current is obtained within the paral- 
lel configuration (positive or negative, depending of the 
specific sign of the polarizations p a ), while for the anti- 
parallel case the spin currents takes less intense values. 
Very interestingly the spin current can become negligible 
also for finite values of Pl/r. Indeed, the above analysis 
of I s shows that when the relation p L « — 3.5p^ is ful- 
filled the spin current goes to zero. In particular when the 
parameters are varied the spin current presents negligible 
values when p L = —kpr, where we have checked that the 
coefficient k depends on the position of energy level e of 
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PL 

FIG. 10: Contourplot of spin current I 3 as a function of the 
polarizations Pl/r by setting the remaining parameters as 

follows: rfr L = 1.5, r£ =J? = i.5, rZ =L/R = 0.3, w = 0.05, 

= 0, yiji = 7r/4, (e — ft) /To ~ 0. Notice that the relation 
J s « Ap L + Bp R with A » 0.94 x 10~ 3 and B ~ 3.3 x 10~ 3 is 
obeyed. 

the QD. Varying it one can rotate the line pi — —npn 
within the (pi,, pn) plane. This property shows an all- 
electrical control of the spin current, potentially useful in 
spintronics devices. 

VI. CONCLUSIONS 

We studied the charge and spin current pumped in an 
interacting quantum dot connected to two external leads 
via tunnel barriers whose transparencies are periodically 
modulated in time via the top gates G1/G2. The cur- 
rent has been calculated by a non-equilibrium Green's 
function (GF) method based on a time dependent slave 
boson approach. In particular, we derived an equation of 
motion of the slave boson field B(t) in terms of an infinite 
series of constraints, generalizing the adiabatic case. By 
using a finite set of constraints (single photon approxi- 
mation) the expressions of the relevant Green's functions 
and of the current have been derived. We showed the 
numerical results of the pumped current for two different 
situations: the spin insensitive and the spin dependent 
case. In the latter case the spin selective properties of 
the system are determined by the presence of two mag- 
netic tunnel barriers. In the spin insensitive case the 
system pumps a charge current that contains beyond a 



term related to the modulation of the out-of-phase ex- 
ternal parameters an additional term due to the internal 
dynamics of the Lagrange multiplier and the slave boson 
field. This term is identified by us as a parasitic pumping 
current. Its presence is due to the phase differences in- 
side the scattering region and is a consequence of strong 
correlations effects. In fact, the charge current, which 
for a free-fermions model depends only on the phase dif- 
ference between the two external parameters, <fn — ipi, 
in the strongly interacting limit depends separately on 
(PL/R. The additional pumping contributions are origi- 
nated by the phase difference of the external parameter 
phase (for instance tpi) and the phase of an internal de- 
gree of freedom (i.e. the slave boson field or the Lagrange 
multiplier) . This mechanism not only works in the inter- 
acting quantum dot case but can operate in the pre sence 
of any internal dynamics, as pointed out in Ref.jlll. \v\. 
Concerning the spin sensitive case the current preserves 
the above characteristics, while the appearance of a non 
vanishing spin current is shown. The spin current is of 
the form I s fts Api + Bpn, where the parameters A, B de- 
pend on the electronic configuration of the QD and pi , pn 
are the polarizations of the magnetic tunnel barriers. By 
considering specific values of the polarization of the bar- 
riers, a change of sign of the spin current can be obtained 
by tuning the energy level e of the QD. This produces a 
rotation of the line I s (pl,Pr) — in the (pl,pr) plane. 
Furthermore, acting on the energy level of the QD the 
system can be driven toward a working point in which the 
pump works as a pure spin current injector. These two 
mechanisms shows the possibility of all-electrical control 
of the spin current pumped through the system relevant 
in spintronics devices. 

Appendix A: Definition of $ (n;,)) (E) 

The definition of the function $(™™(.E) in terms of the 
GFs of the quantum dot region takes the following form: 

= £ / J^WW (Al) 

x [f(E 2 -nu})G r a . rl (E + E 2 -nu) 
+ ^G< r ,(E + E 2 -ncj)}, 

where the function is defined as follows: 

= r a , B (A2) 

^=±1) = r CT;0 B ± i + r CT;±1 So, 

With r CT:n = Ea r ";n- 
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